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ABSTRACT 

Using the most recent observational data concerning the Extragalactic Back- 
ground Light and the Radio Background, for a source at an arbitrary redshift in the 
range 1CP 6 < z s < 3 we compute the energy Eq of an observed 7-ray photon in the 
range 10 GeV < E a < 10 13 GeV such that the resulting optical depth t^(E , z s ) takes 
the values 1, 2, 3 and 4.6, corresponding to an observed flux dimming of e _1 ~ 0.37, 
e~ 2 ~ 0.14, e~ 3 ~ 0.05 and e~ 4 - 6 ~ 0.01, respectively. Below z s ~ 10~ 6 we find that 
t 7 (Eq,z s ) < 1 for any value of E . In the limiting case of a local Universe (z s ~ 0) 
we compare our result with the one derived in 1997 by Coppi and Aharonian. The 
present achievement is of paramount relevance for the planned ground-based detectors 
like CTA, HAWC and HiSCORE. 

Key words: diffuse radiation - gamma-rays: observations. 



1 INTRODUCTION 

Very-high-energy (VHE) astrophysics is on the verge to en- 
ter its golden age. Planned ground-based detectors like CTA 
(Cherenkov Telescope Array) (Actis et al. 2011), HAWC 
(High Altitude Water Cherenkov Experiment) (Sinnis 2005) 
and HiSCORE (Hundred Square-km Cosmic ORigin Ex- 
plorer) (Tluczykont et al. 2011) will probe within the next 
few years the energy range from 10 GeV up to 10 GeV (CTA 
and HAWC) and even up to 10 9 GeV (HiSCORE) with un- 
precedented sensitivity. 

Unfortunately, a stumbling block along this exciting av- 
enue is the existence of a soft photon background in the Uni- 
verse which leads to a strong suppression of the observed flux 
through the 77 — > e + e~ pair-production process. The onset 
of this process depends both on the observed energy and on 
the source redshift, as it will become apparent later. Specifi- 
cally, it is mostly due to the Extragalactic Background Light 
(EBL) - i.e. to the background in the infrared, visible and 
ultraviolet region - in the energy range 10 GeV — 10 5 GeV, 
to the Cosmic Microwave Background (CMB) in the range 
10 5 GeV - 10 10 GeV and to the Radio Background (RB) 
in the range 10 10 GeV - 10 13 GeV. Moreover, even though 
the spectral energy distribution (SED) of the EBL has re- 
mained quite uncertain for a long time, a remarkable agree- 
ment among the various EBL models has been reached re- 
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cently, and also the RB has been measured in 2008 with 
considerably better precision than before. 

All this prompts us to carefully evaluate the optical 
depth - which quantifies the photon absorption due to 
the above pair-production process - in the energy range 
10 GeV — 10 13 GeV for a source redshift up to z s = 3. As 
a particular case, we derive the photon mean free path for 
77 — > e + e~ as a function of energy in the local Universe 
(z s ^ 0) in order to compare it with that obtained by Coppi 
and Aharonian (CA) in 1997 (Coppi & Aharonian 1997). A 
less detailed but similar result was independently derived at 
the same time by Phrotheroe and Biermann (Phrotheroe & 
Biermann 1997). 



2 EVALUATION OF THE PAIR-PRODUCTION 
CROSS-SECTION 

We start by recalling that the photon survival probability 
P 1 ^ 1 (Eo, z s ) is currently parameterized as 

P^(Eo,z a ) = e - T ^ E °' z ° ) , (1) 

where Eo is the observed energy and t 7 (-Eo,2 s ) is the opti- 
cal depth that quantifies the dimming of the source at red- 
shift z s . Clearly r 7 (i?o, z s ) increases with z s , since a greater 
source distance implies a larger probability for a photon 
to disappear. Apart from atmospheric effects, one typically 
has Ty(Eo,z a ) < 1 for z s not too large, in which case the 
Universe is optically thin all the way out to the source. 



2 De Angelis, Galanti, Roncadelli 



But depending on Eo it can happen that t^(Eo,z s ) > 1, 
and so at some point the Universe becomes optically thick 
along the line of sight to the source. The value Zh such that 
t 7 (15o, Zh) = 1 defines the "/-ray horizon for a given Eo, and 
sources beyond the horizon (namely with z 3 > Zh) tend to 
become progressively invisible as z 3 further increases past 

Zh- 

Whenever dust effects can be neglected, photon deple- 
tion arises solely when hard photons of energy E scatter off 
soft background photons of energy e permeating the Uni- 
verse, which gives rise to hard photon absorption. Let us 
proceed to quantify this issue. 

Regarding E as an independent variable, the process is 
kinematically allowed for 



e > e t hr(E, ip) 



(2) 



E (1 — cos ip) 

where ip denotes the scattering angle and m e is the electron 
mass. Note that E and e change along the line of sight in 
proportion of 1+2 because of the cosmic expansion. The cor- 
responding Breit- Wheeler cross-section is (Breit & Wheeler 
1934; Heitler 1960) 
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with 



W(0) = (1 - p 2 ) 
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The cross-section depends on E, e and ip only through the 
speed j8 — in natural units - of the electron and of the 
positron in the center-of-mass 

r, 2 4 -|V2 



P(E,e,tp) = 



1 - 



Ee (1 — cos ip) 



(4) 



and Eq. S implies that the process is kinematically al- 
lowed for /3 2 > 0. The cross-section a- 11 {E,e,p) reaches its 
maximum cr™ 7 ax ~ 1.70 • 10" 25 cm 2 for /3 ~ 0.70. Assuming 
head-on collisions for definiteness (ip — n), it follows that 
a 11 {E,e,ii) gets maximized for the background photon en- 
ergy 



cV 



(5) 



where E and e correspond to the same redshift. Explicitly, 
the situation can be summarized as follows. 

• For 10 GeV < E < 10 5 GeV the EBL plays the leading 
role. In particular, for E ~ 10 GeV a 11 (E, e,n) is maximal 
for e ~ 50 eV corresponding to far-ultraviolet soft photons, 
whereas for E ~ 10 GeV a 11 (E, e, ir) is maximal for e ~ 
5- 10~ 3 eV corresponding to soft photons in the far-infrared. 

• For 10 5 GeV < E < 10 10 GeV the interaction with the 
CMB becomes dominant. 

• For E > 10 10 GeV the main source of opacity of the 
Universe is the RB. 



3 EVALUATION OF THE OPTICAL DEPTH 

Within the standard ACDM cosmological model r 7 (_Eo,2„) 
arises by first convolving the spectral number density 



n 1 (e(z), z) of background photons at a generic redshift z 
with a 11 (E(z), e(z), ip) for fixed values of z, tp and e(z), and 
next integrating over all these variables (Gould & Schreder 
1967; Fazio & Stecker 1970). Hence, we have 



y(E , 



dz 
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d(cos tp) 
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: / de(z)n-,(e(z),z)o-^(E(z),e(z),tp) , 

where the distance travelled by a photon per unit redshift 
at redshift z is given by 



dl(z) 
dz 



H o (1 + z) [fi A + flu (1 + zf] 
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(7) 



with Hubble constant Hq ~ 70 km s _1 Mpc -1 , while J1a — 
0.7 and Qm — 0.3 represent the average cosmic density of 
matter and dark energy, respectively, in units of the critical 
density p cr ~ 0.97- 10" 29 gcm- 3 . 

Once n~f(e(z), z) is known, t 7 (_Eo,z) can be computed 
exactly; generally the integration over e(z) in Eq. Q must 
be performed numerically. 

Finally, in order to get a feeling about the considered 
physical situation, it looks suitable to discard cosmologi- 
cal effects, which evidently makes sense only for z s small 
enough. Accordingly, z 3 is best expressed in terms of the 
source distance D — cz s / Hq and the optical depth becomes^ 

where A 7 (B) is the photon mean free path for 77 — > e + e~ 
referring to the present cosmic epoch. As a consequence, Eq. 
|l| becomes 

P^(E,D) = e - D/x ^ E) . (9) 



4 SOFT PHOTON BACKGROUND 

Our main goal is at this point the determination of 
n-y(e(z),z). For the sake of clarity, we consider separately 
the EBL, the CMB and the RB. 

The EBL density n 7 (e(z),z) is in principle affected by 
large uncertainties arising mainly from foreground contami- 
nation produced by zodiacal light which is various orders of 
magnitude larger than the EBL itself (Hauser & Dwek 2001). 
Below, we sketch schematically the different approaches that 
have been pursued, without any pretension of completeness. 

• Forward evolution - This is the most ambitious ap- 
proach, since it starts from first principles, namely from 
semi-analytic models of galaxy formation in order to predict 
the time evolution of the galaxy luminosity function (Pri- 
mack et al. 2001; Primack, Bullock & Sommerville 2005; 
Gilmore et al. 2009; Gilmore et al. 2012). 

• Backward evolution - This begins from observations 
of the present galaxy population and extrapolates the 
galaxy luminosity function backward in time. Among oth- 
ers, this strategy has been followed by Stecker, Malkan 
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and Scully (Stecker, Malkan & Scully, 2006) - whose re- 
sult has unfortunately been ruled out by the measurements 
by Fermi/LAT (Abdo et al. 2010) - and by Franceschini, 
Rodighiero and Vaccari (FRV) (Franceschini, Rodighiero & 
Vaccari 2008). 

• Inferred evolution - This approach models the EBL by 
using quantities like the star formation rate, the initial mass 
function and the dust extinction as inferred from observa- 
tions (Kneiske, Mannheim & Hartmann 2002; Kneiske et al. 
2004; Finke, Razzaque & Dermer 2010). 

• Observed evolution - This method developed by 
Dominguez and collaborators (D) relies on observations by 
using a very rich sample of galaxies extending over the red- 
shift range < z s < 2 (Dominguez et al. 2011). 

• Compared observations - This technique has been im- 
plemented in two different ways. One consists in compar- 
ing observations of the EBL itself with blazar observations 
with Imaging Atmospheric Cherenkov Telescopes (IACTs) 
and deducing the EBL level from the VHE photon dim- 
ming (Schrodter 2005; Aharonian et al. 2006a; Mazin & 
Raue 2007; Mazin & Gobel 2007; Finke & Razzaque 2009). 
The other starts from some 7-ray observations of a given 
blazar below 100 GeV where EBL absorption is negligible - 
typically using Fermi/LAT data - and infers the EBL level 
by comparing the IACT observations of the same blazar with 
the source spectrum as extrapolated from former observa- 
tions (Orr, Krennrich & Dwek 2011) (but see also Costa- 
mante 2012). In the latter case the main assumption is that 
the emission mechanism is presumed to be determined with 
great accuracy. In either case, the crucial unstated assump- 
tion is that photon propagation in the VHE band is governed 
by conventional physics. 

• Minimal EBL model - Its aim is to provide a strict 
lower limit on the EBL level. It relies on the same strategy 
underlying the inferred evolution, but with the parameters 
tuned in such a way to reproduce the EBL measurements 
from galaxy counts (Kneiske & Dole 2010). 

Quite remarkably, all methods - apart obviously the last 
one - yield basically the same results in the redshift range 
where they overlap, so that at variance with the time when 
the CA analysis was done (1997) nowadays the SED of the 
EBL is fixed to a very good extent. We will employ here the 
FRV model (Franceschini, Rodighiero & Vaccari 2008) but 
we shall check our results using the D model (Dominguez 
et al. 2011). As far as the CMB is concerned we take the 
standard temperature value T — 2.73 K. Finally, the most 
recent available data for the RB are employed (Gervasi et 
al. 2008), with a low-frequency cutoff taken at 2 Mhz. 



5 RESULTS 

Proceeding according to the above strategy - namely tak- 
ing into account the EBL, the CMB and the RB - we di- 
rectly evaluate T 7 (i5o,Zs) over the energy range 10 GeV < 
E < 10 13 GeV and within the source redshift interval 
10 -3 < z a < 3, which is the range where the FRV model 
yields the contribution of the EBL to r 7 (_Eo,z 3 ) (the D 
model is restricted to 10~ 2 < z s < 2). In the lack of any in- 
formation about the EBL for z s < 10~ 3 , we linearly extrap- 
olate Tj(Eo,z s ) down to z s — 10~ 6 . Such an extrapolation 
looks quite reliable not only since r 7 (_Bo, z 3 ) behaves linearly 
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Figure 1. Source redshifts z s at which the optical depth takes 
fixed values as a function of the observed hard photon energy Eq . 
The curves from bottom to top correspond to a photon survival 
probability of e _1 ~ 0.37 (the horizon), e~ 2 ~ 0.14, e -3 ~ 0.05 
and e~ 4,6 ~ 0.01. For z s < 10 -6 the photon survival probability 
is larger than 0.37 for any value of Eq. 



2 but mainly because 
and so Eq. (JsJ indeed 



already in the range 10 < z s < 10 
at such low redshifts we have D oc z s 

implies r 7 (_Eo, z s ) oc z s . Moreover, we stress that the extrap- 
olation in question does not practically affect our result. For, 
as it is evident from Fig. [I for z s < 10 -3 the energies Eq 
for which Tj(Eo, z s ) takes our prescribed values tend to ex- 
ceed 10 5 GeV, which means that our result for z s < 10~ 3 
is dominated by the CMB and the RB rather than by the 
EBL. Although we have computed r 7 (J5o, z 3 ) using the FRV 
model, we have checked that it basically remains unaffected 
by employing the D model in the redshift range where they 
overlap. This is our main result, which is plotted in Fig. 
[I] where the solid line corresponds to T~,(Eo,z a ) = 1, the 
dot-dashed line corresponds to r 7 (_Eo,z s ) = 2, the dashed 
line corresponds to r 7 (_Eo,2 s ) = 3 and the dotted line cor- 
responds to r 7 (_Eo, z a ) = 4.6, which give rise to an observed 
flux dimming of about 0.37, 0.14, 0.05 and 0.01, respectively. 
We see that for z s < 10 -6 it turns out that Tj(Eo,z s ) < 1 
for any value of Eq . 

In order to compare our finding with the CA result, we 
disregard cosmological effects thereby computing the 77 — > 
e + e~ mean free path A 7 (_E) in the local Universe (z s — 0) 
by using Eq. Q with Tj(Eq, z s ) evaluated by means of Eq. 
^ with z s = 10" 6 . 

To facilitate the comparison between our and the CA 
results we have superposed in Fig. [5] the behavior of A 7 (E) 
as found with the above procedure - represented by red 
dot-dashed line - over Fig. 1 of CA, where the black lines 
represent the similar results derived by them (see captions). 
There is of course little surprise that at z s — our result 
is barely in agreement with that of CA at z s = 0. Still, 
it should be appreciated that the improved behaviour of 
A 7 (-B) found here is in disagreement with the one arising 
from a single choice of the EBL model in the CA analysis. In 
addition, while the CMB contribution is obviously identical 
in both cases, as far as the RB is concerned our curve - 
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Figure 2. Comparison of the mean free path A 7 for 77 — > e+e _ 
as derived in the text using a low-frequency cutoff at 2 Mhz (red 
dot-dashed line) with the one obtained by CA in the local Uni- 
verse (z a = 0) as a function of the observed energy E. The CA 
black lines labelled by a, b, c represent three different EBL models, 
the black solid line corresponds to the CMB, the black lines la- 
belled by 1, 2, 3 represent a model of the RB with low-frequency 
cutoff at 5, 2 and 1 Mhz, respectively, and the black triangles 
correspond to the RB under the assumption that it is completely 
extragalactic. 

which corresponds to a low-frequency cutoff of 2 Mhz - lies 
between those corresponding to the low-frequency cutoffs of 
2 Mhz and 1 Mhz, respectively, of the CA result. 



6 CONCLUSIONS 

We have quantified by means of the optical depth the pho- 
ton absorption caused by the pair-production process in the 
observed energy range 10 GeV < E < 10 13 GeV and for a 
source redshift up to z s = 3. We have found that depending 
on Eo the absorption can be quite large for z B > 2 ■ 10 -6 and 
becomes dramatic around 10 6 GeV. However, for z B < 10~ 6 
the absorption is irrelevant for any value of Eo . 

As it is clear from Fig. [TJ our conclusion is of great 
importance for the planned VHE detectors like CTA, HAWC 
and HiSCORE. 

Moreover, we have specialized our analysis to the local 
Universe (z s ~ 0) where the optical depth is more conve- 
niently replaced by the mean free path, and we have com- 
pared it with the same quantity as evaluated by Coppi and 
Aharonian in 1997. 

It has been realized that the blazars observed so far 
by IACTs give rise to the pair-production anomaly. The 
H.E.S.S. collaboration (Aharonian et al. 2006b) first ob- 
served that the SED from the blazars H 2356-309 and 1ES 



1101-232, at redshifts z s — 0.165 and z s = 0.186, respec- 
tively, could be explained by a very low EBL level; subse- 
quently, the MAGIC collaboration (MAGIC Collaboration 
2008) have reinforced this evidence with the data from the 
AGN 3C279 at z B = 0.54. Later, De Angelis et al. (De An- 
gelis et al. 2009) have observed this effect in the spectral 
indices at VHE of a sample of AGN at z a > 0.2. Recently a 
statistical analysis of the SED from all blazars observed at 
VHE indicates a level of the EBL lower than that predicted 
even by the minimal EBL model; the indication remains at 
a confidence level between 2.6 and 4.3 depending on the 
adopted EBL model (Horns & Meyer 2012; Meyer, Horns & 
Raue 2012). Among many other things, we find it very in- 
teresting to see whether this effect persists also at energies 
much higher than those presently explored by IACTs, and 
our result looks essential as a benchmark for comparison. 

In a forthcoming paper an analysis along the same lines 
of this work will be carried out taking Axion-like particles 
(ALPs) into account, since photon- ALP oscillations tend to 
drastically reduce photon absorption effects of the kind con- 
sidered here, thereby considerably enlarging the 7-ray hori- 
zon at the VHEs that CTA, HAWC and HiSCORE will be 
able to probe (see De Angelis, Galanti & Roncadelli 2011 
and references therein). Moreover, it has very recently been 
pointed out that such a mechanism would explain the pair- 
production anomaly in a natural fashion and works for val- 
ues of the photon-ALP coupling in the reach of the planned 
upgrade of the ALPS experiment at DESY (Meyer, Horns 
& Raue 2013). 
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